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Abstract. A complete solution to the multiplier version of the inverse problem of the calcu- 
lus of variations is given for a class of hyperbolic systems of second-order partial differential 
equations in two independent variables. The necessary and sufficient algebraic and differential 
conditions for the existence of a variational multiplier are derived, ft is shown that the number of 
independent variational multipliers is determined by the nullity of a completely algebraic system 
of equations associated to the given system of partial differential equations. An algorithm for 
solving the inverse problem is demonstrated on several examples. Systems of second-order par- 
tial differential equations in two independent and dependent variables are studied and systems 
which have more than one variational formulation are classified up to contact equivalence. 



1. Introduction 

In this paper we examine the inverse problem of the calculus of variations for systems of 
partial differential equations 

(1.1) u% y = f a (x,y,u 7 ,u2,u] / ), a, 7 = 1,2,..., m, 

in m dependent variables u a and two independent variables x and y. Systems (jl.ip arise in 
various contexts such as symmetry reduction in general relativity non-linear cr-models, and 
generalized Toda lattices (see [11] , [13] , [20] , [21] , and [23] ) . For systems (jl.ip we give a complete 
solution to the variational multiplier version of the inverse problem. In particular, we give an 
explicit algorithm for determining the number of possible inequivalent Lagrangians for a given 
system (jl.ip . 

The variational multiplier problem is stated as follows: Given a system of differential equations 
F a (x' 1 , u 7 , u[, u]j,ujj k , . . . ) =0, does there exist a Lagrangian L(x l , u 1 ,u] , ujj, . . . ) and functions 
M%, with det (Ml) ^ 0, such that 

(1.2) E a (L) = M%F 1 , 

where E is the Euler-Lagrange operator. If such an M" exists, then it is referred to as a 
variational multiplier. If (|1.2j) is satisfied, the equations F a = are equivalent to a system of 
Euler-Lagrange equations in the sense that solutions of the Euler-Lagrange equations for L are 
solutions to F a = and conversely solutions to F a = are solutions for E a {L) = 0. 

Also of importance is the number of Lagrangians for a particular system. It is well known 
that two Lagrangians L and L have identical Euler-Lagrange expressions if and only if L and 
L differ by a total divergence, at least locally. Therefore we consider two Lagrangians to be 
equivalent if they differ by a total divergence. We note that it is possible to have a system of 
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differential equations which have two or more inequivalent Lagrangians (and thus more than one 
multiplier). 

The variational multiplier problem is of some interest in theoretical physics. It is widely ac- 
cepted that fundamental physical theories can be derived from an action principle, or Lagrangian. 
For a given theory it is important to determine whether the action is unique. Examples with 
multiple Lagrangians exist in Newtonian Mechanics and the SU(2) chiral model ([HJ, [15),) 

The simplest case of the variational multiplier problem is for a scalar second-order ordinary 
differential equation u xx = F(x,u,u x ). It has been shown by several authors, including Dar- 
boux [10], that any scalar second-order ODE is variational and the most general multiplier 
(and Lagrangian) depends on an arbitrary function of two variables. The multiplier problem 
for systems of second-order ordinary differential equations has been studied by many authors 
including Douglas [9], Anderson and Thompson [3], Thompson, Crimpin, Sarlet, Prince and 
Martinez (See [7] , [8] , [22] , and [23].) Our solution to the multiplier problem for systems (jl.ip 
is partially based on ideas developed in Anderson and Thompson's paper [3J. In that paper 
Anderson and Thompson used the variational bicomplex to derive a system of algebraic and 
differential conditions, with components of the multiplier matrix as unknowns, for the existence 
of a multiplier. They showed that there is a one-to-one correspondence between variational 
multipliers for the given system and certain cohomology classes in the variational bicomplex 
associated to the given system of differential equations. While a significant amount of research 
has been done on the variational multiplier problem for second-order ODE systems, a complete 
solution remains elusive in the sense that there is no general closed form characterization, in 
terms of invariants for the system, for determining the existence and degree of uniqueness of 
multipliers for the given system. 

The variational multiplier problem for higher order scalar ordinary differential equations 
has been studied by Fels [12] and Juras [17J. Fels [12] obtained a complete solution to the 
variational multiplier problem for fourth-order scalar ordinary differential equations u xxxx = 
F(x ,u,u x ,u xx ,u xxx ). Using Cartan's method of equivalence, he was able to produce two differ- 
ential invariants whose vanishing completely characterizes the existence of a variational multi- 
plier. Unlike the second-order case, the multiplier is unique up to a constant multiple. In |17| . 
Juras obtained a similar solution for sixth and eighth-order scalar ordinary differential equations, 
although the differential invariants are increasingly complicated for higher order systems. 

For partial differential equations, there are fewer papers on the variational multiplier problem. 
Anderson and Duchamp [2] studied the variational multiplier problem for scalar second-order 
quasilinear partial differential equations F = where F = (x k ,u,Uk)uij + B(x k , u, Uk)- 
Anderson and Duchamp [2] proved that if det(A l:) ) ^ 0, then F has a variational multiplier if 
and only if a certain 1-form \ ls closed. Moreover, x 1S expressed explicitly in terms of F and 
its derivatives and the multiplier, if it exists, is unique up to a constant multiple. They also 
showed that second-order scalar evolution equations are never variational. Juras [16] examined 
the inverse problem for a scalar hyperbolic second-order partial differential equations in two 
independent variables. He was able to show that an equation is variational if and only if two 
particular differential invariants H and K are identically equal. Moreover, the multiplier is 
unique up to a constant. Juras' result is equivalent to the Anderson-Duchamp result if the 
equation is quasilinear. However, his results also apply to hyperbolic Monge- Ampere equations. 

In this paper we study the variational multiplier problem for /-Gordon systems 

(1.3) u^y = f a {x,y,u i ;ul,u'l), a, 7 = 1, 2, . . . , m. 

We remark that we are working under the assumption that the functions f a are C°° on some 
open set U C R 2 x R 3m . In the following section we outline the complete solution to the 
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multiplier version of the inverse problem for systems fjl .3[) . In particular, we give an algorithm 
for determining the number of variation multipliers (and Lagrangians) for a given system (|l,3p . 
We show that the most general multiplier depends on finitely many constants determined by 
the dimension of the nullspace of a certain matrix depending on the functions f a and their 
derivatives. In Section 3 we demonstrate our algorithm for solving the inverse problem on several 
examples. In Section 4 we classify all systems (jl.3p in two dependent variables that admit two 
or more inequivalent Lagrangians. In Sections 5 and 6, we turn our attention to proving two 
technical propositions stated in Section 2. In particular, we define the variational bicomplex for 
systems (|1 .3j) and show there is a one-to-one correspondence between the Lagrangians for systems 
(|1.3p and special classes of 3-forms u with duj = 0. We then derive necessary and sufficient 
algebraic and differential conditions for the existence of a variational multiplier. In the appendix 
we prove a general theorem on the existence of solutions to systems of combined differential and 
algebraic equations that allows us to determine the number of variational multipliers from purely 
algebraic data. 

This work is an extension of a result established in my PhD dissertation. I wish to thank 
my advisor Ian Anderson for suggesting this problem and an uncountable number of helpful 
discussions on this subject. 



2. Main Results 

In this section we state our solution to the variational multiplier problem for systems (jl.lj) . 
Our solution depends on two propositions which we will prove in Section 6. The first proposition 
gives us a general normal form for variational systems (jl.ip and states that the Lagrangian 
associated to a particular variational multiplier is unique up to modification by a total divergence. 

Proposition 2.1. If there exists a first- order Lagrangian L(x, y, u 7 , u%, Uy) and a non- degenerate 
variational multiplier M aj g(:c,y, u 7 ,uJ,Uy) such that 

E a {L) = M ap (<, - f(x, y, u\ u%,ul)) , 

then M ap = M a p(x, y, u 7 ), M aj3 = M 0a , 

q2 fa Q2 fa 

(2.1) — J—=0 — ^_^ = o, Va,/3,7 = 1,2,... ,m. 

OUxOUx OUyOUy 

Moreover, if L(x,y,u Y ,uZ,Uy) is another Lagrangian associated to the variational multiplier M a @ 
such that E a (L) = M a/3 {vP xy — f 13 ), then L = L + Div Q, where Q = ( Qi(x, y, u 7 ), Q2(x, y, u<) ) 
and Div Q = D X Q\ + D y Q2- 



It follows immediately from Proposition (|2.ip that we may restrict ourselves to systems of the 
form 

(2.2) u% + C^(x, y, u T ) ui n 7 + A%(x, y, u> 7 + B«(x, y, u> 7 + E a (x, y, u T ) = 0. 

We remark that the condition that M a p is symmetric and has no first-order derivative depen- 
dence also follows from a result established by Henneaux in [14J. 
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Of paramount importance in our study of the inverse problem are three quantities HZ, KZ, 

a/3 



and S^r defined by 



(23a ) jt - dP i dP dF d2f " d2p -J 92 P f Q2P -J 

{ ' a du« + du°du« du-dx du%dvP x du«d4 J du^dui xx ' 

(23b) W= dP I dP dP d2f " d2p -J d2p f g 92P -J 

1 ' ; a dug du« du«d y du«duP y du?du p J du a du p yy ■ 



(2.3c) & 



d 2 p d 2 p 



•a 

01 ~ du£du° du%du§ 



In the formulas (|2.3|) and throughout this paper we are adopting the Einstein summation con- 
vention where repeated indices are summed upon. We remark that the quantities HZ, KZ, and 
(Slg are relative invariants under the pseudo-group of local contact transformations that pre- 
serve systems (jl.ip . The functions H" and are natural extensions of the generalized Laplace 
invariants defined in [3] for scalar hyperbolic equations. For systems (|2.2p it follows from (|2.3p 



that the functions HZ and KZ have no second-order derivative dependence and SZs has no first- 
order derivative dependence, that is HZ = HZ(x,y,u T ,u%,Uy), KZ = KZ[x,y,u T ,u^,Uy), and 
S a = SZf3{x,y,u T ), where HZ and KZ are quadratic in the variables u@ and u@. Our second 
proposition characterizes the algebraic and differential conditions on a variational multiplier for 
a system ()2.2|) . 

Proposition 2.2. A system (12 - 2 j) is multiplier variational with a multiplier M a p if and only if 
M a/3 satisfies det M aB ^ and 

(2.4a) M aa Hp = MftoKZ, M aa S^ = -M^S^, 

(2.4b) dM ap = M aa % + M &a 122, 

where fi£ = C a ar du T +A£dy + B a a dx. 

Remark 2.3. In Section 6, we will derive the algebraic and differential conditions (|2,4p and 
give an algorithm for constructing the first-order Lagrangian L associated to a variational mul- 
tiplier M a p satisfying (|2.4p . The algebraic conditions (|2.4ap include the integrability conditions 
for d 2 M a(3 = for fObl) and genuine algebraic constraints not arising from the integrability 
conditions. We will show that the set of all solutions to (|2.4p is a finite dimensional real vector 
space with the dimension determined by the nullity of a completely algebraic system of equa- 
tions. In particular, for any system of the form (|2.2p . the solutions to (|2.4p are completely 
determined. However, in some cases there are non-trivial solutions M a/3 to (12. 4p that do not 
satisfy det M a p ^ 0. In practice, we first determine a basis for the solutions to (12.4p and then de- 
termine if a non-degenerate variational multiplier can be constructed from a linear combination 
of the basis solutions. 

We are now ready to state our solution to the inverse problem which says that the variational 
multipliers M af3 satisfying (12. 4p are completely characterized by the solutions to an algebraic 
system $"' 3 M a p = 0, where <3? depends on the given system u" y = f a . 

Theorem 2.4. Let uf. y = f a (x,y,u l ,ul,Uy) be a system of m partial differential equations 
where f a £ C°°(U) for some open set U C R 2 x R 3m . Then there is a matrix 3>, with Rank(&) < 
m(m + l)/2, depending the functions HZ, KZ, and S 1 ^ and their derivatives, whose nullspace 
completely determines the number of linearly independent solutions to (|2.4p . Specifically, if 
r = Rank 5>(z) is constant for all points z in a open set V C U, then at every point zq G V 
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there is a neighborhood W C V of zq where the set of solutions to (j2.4[) is a m(m + l)/2 — r 
dimensional vector space over R. 

Proof. The algebraic system for a multiplier M af3 is constructed as follows. As a consequence 
of proposition (|2,ip . the functions M a p have no 1-jet dependence and Ha and K2 are quadratic 
in u € x and u r y . We can decompose the algebraic condition M a ^Hj = M^K2 into several linear 
algebraic systems on the components of the multiplier M af3 , with the coefficients depending only 
on x, y, and u. We then express all algebraic conditions (|2.4ap on M Q7 as a single system of 
linear equations 

(2.5) ($ )f (x,y,u T )M a0 = 0, a =1,2,...,**,, 

where may be viewed as fco x m(m + l)/2 matrix. Differentiating (|2.5p and substituting from 
(|2.4bp . we get the algebraic condition 

(2-6) [ + (*o)r ^ + ( $ o)f O?] M a/3 = 0. 

Equation (|2.6[) represents a system of purely algebraic conditions on M af3 . We then express the 
combined algebraic systems ()2.5p and (|2.6p as 



(*i)fM a/J = J o = l,...,fci. 
We can proceed inductively to define a system of equations 

(2.7) M a/3 = 0, a= l,...,fc i+ i, 
where (I2.7P consists of the system (^J^M^ = along with the system 

[(d$i)f + 2(^)r^]M a ^ = 0. 

Clearly, at each point zo = (^0)2/0) ^q) £ U C R 2 x R m and for all j > 0, the ranks of the 
matrices <3?j satisfy 

< Rank $0(^0) < Rank *i(«o) < • • • < R ank $j( z o) < — ^ -• 

If given a point zq € [7, we see that after k < m(m + l)/2 differentiations, the rank of the 
matrices $>i(zo) must stabilize at some < I < m(m + l)/2. More precisely, at the point zq we 
have 

I = Rank $k( z o) = Rank <&k +i (zo), Vi > 0. 

Let {C^ 7 , , . . . , C^, 7 }, where s = m(rn + l)/2 — I, be a basis for the set of solutions to the 
system of linear equations 

(2.8) ($ fc )^M a7 = 

at the point zq. In particular, {C^, 7 | i = 1, . . . } is a linearly independent collection of constant 
bilinear forms such that k)""' (zo) C ai = 0. 

If the rank of is constant in a neighborhood V C Uofzg, then a general result on systems of 
algebraic-differential equations, which we prove in the appendix, states that there exists a neigh- 
borhood W C V of zq and a linearly independent collection of functions {M^p, M^p, ■ ■ ■ , -M^} 
such that 

(2.9) M^(z o ) = C i a0 , ($ k )f(z)M a0 (z) = O, dM ap = M aa n° + M^U*, 
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for all z £ W, and i = 1, 2, . . . , s. Moreover, we claim that if a given collection of functions M a/3 
satisfies the algebraic and differential conditions (j2.9[) . then M a p can be expressed as a linear 
combination 

s 

(2.10) M ap = Y,CiM l af3 , Ci €R. 

i=l 

Indeed, if M af3 satisfies the algebraic condition (^k)^ (z)M a p{z) = for all z € W, it follows 
that M Q(3 can be uniquely expressed as 

(2.11) M a0 (z) = J2a(z)K f3 (z). 

If M a/3 satisfies the differential condition (|2.4bl) . then it follows from (12. 9j) and (|2.1ip that 
M % a p dci = 0. Since the functions MIg are pointwise linearly independent in some neighborhood 
of 2q, we deduce that dci = 0, which in turn implies that a € R for z = 1,2,... , s. This 
establishes fj2. 10f> and completes the proof of the theorem. □ 



3. Examples 



In this section we demonstrate our algorithm for solving the variational multiplier problem 
on several examples. For each example we calculate the invariants and using 

(|2.3h and we explicitly list the initial algebraic conditions H2.4aH on a multiplier M a p. We then 
determine the differential condition (|2,4bl) and differentiate the algebraic conditions to uncover 
any additional algebraic constraints. In each case we find the most general Lagrangian and 
multiplier for the given system. 



Example 3.1. For our first example, consider the system 
(3-1) 



u. 



^xy — u : u xy 

We will show that (|3.ip admits two Lagrangians. In this case 
algebraic condition from (|2.4p is 



0, so that the only nontrivial 



(3.2) 

The differential condition is dM a t 



0. Using (|2.3p we calculate 



H = K 



1 

1 



and it follows from (|3.2p that the only algebraic constraint is M\\ = M-2i- Since the differential 
condition is dM af3 = 0, differentiating Mu = M22 produces no additional algebraic conditions. 
The most general multiplier in this case is 



M 



a b 
b a 



a, b 6 R, a 2 - b 2 ^ 0. 



The Lagrangian corresponding to the multiplier M is given by 



L 



' : (U X Uy + V X Vy + 2UV) 



-(2u x Vy + u 2 + v 2 ) 



A routine calculation shows that the Euler-Lagrange equations for L are 



E\{L) = a(u xy — v) + b(v, 



J xy 



E 2 {L) = b(u xy - v) + a(v. 



Jxy 



U 
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Example 3.2. Consider the system 

(3.3) u xy = v, v xy = xu. 

The only difference between this system and the one in the first example is the x in the second 
equation. We will show that (|3.3|) admits a unique Lagrangian. The initial algebraic conditions 
(|2.4p on the multiplier M a/3 reduce to 

(3.4) M n = xM 22 . 

The differential condition is again dM a/ 3 = 0. Differentiating (13, 4p and substituting dM a p = 
results in M 22 dx = 0, implying that M22 = Mn = 0. Further differentiations add no new 
algebraic conditions. Consequently the multiplier M, where Mn = M 22 = and M\ 2 = M 2 \ = 
1, is unique up to scalar multiplication. The Lagrangian for (|3.3[) is L = —u x v y — (u 2 + xv 2 )/2. 

Example 3.3. For our third example, we again we make a slight change on the first example 
(|3.1|) to get a system that is not variational. Let 

(3.5) Uxy = V, V X y = U x . 

In this case the algebraic and differential conditions (|2.4p on the multiplier M Q( g are given by 

(3.6) Mn = 0, dM n = 2M 12 dy, dM 12 = M 22 dy, dM 22 = 0. 

Differentiating the algebraic condition Mn = we find that M12 dy = 0, implying that M12 = 0. 
Differentiating M12 = 0, it follows that M22 dy = 0. Consequently, the only solution to (|3,6p is 
the trivial solution M = 0. 

Example 3.4. Consider the system 

(3.7) M ; + r; 7 ( M e )^ = o 1 

where r^* are the components of a symmetric connection on an m-dimensional manifold Ai. We 
will show that (|3.7p is variational if and only if T is a metric connection. Since T is symmetric 
we have 5^ = and the Laplace invariants for (|3.7p are given by H2 = R e J a aU%Uy and 
Ka = Ra aeU^Uy, where R e a 1/3 are the components of the curvature tensor associated to the 
connection T. Using properties of the curvature tensor, we can show that the conditions (|2,4p 
on the multiplier are 

(3.8) M ai R p \ a + M^R a \ a = 0, dM ap = (M ai T p \ + M^T^ t ) du e . 

We see immediately from the differential condition in (|3.8p that any multiplier M a/3 must sat- 
isfy dMap/dx = and dM a p/dy = 0. Consequently, M a p = M a p{u e ). It follows that the 
differential condition on the multiplier simplifies to V 7 M aj a = 0, where V 7 denotes covariant 
differentiation with respect to u 1 . This proves that T is the Levi-Cevita connection for the metric 
M a p. Subsequent differentiations of the algebraic condition (|3.8p imply that M a/3 must satisfy 
M a7 Vj -R 7 e/3o - = M(3 7 V/ R?aae, where V/ = V Ql V Q2 . . . V Qfc . We remark that if (M, T) is a (lo- 
cally) symmetric space, then the algebraic condition (|3.8p involving the curvature completely 
determines the number of linearly independent metrics for the given connection. 

Example 3.5. Consider the system of differential equations 

(3.9) u% + c a p y x ui = ^ 

where C" 7 G R are the structure constants of an m dimensional Lie algebra 0. The constants 
Cg are skew symmetric in the lower indices and satisfy the Jacobi identity. Using (|2.3|) we 
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calculate the invariants 



"a 



Q7 
°a0 



x^y i 



C^cfi 'er) U x U y 1 



°/3o 



Consequently, the algebraic conditions (12.-4H can be summarized as 



M Q7 (C2 a Cp, 



0T 



/"(7 /HKT \ 

i7 



M 07 (CJ T Cf Q 



(3.10a) 

(3.10b) M ai C} e + Af^C^ = 0. 

The differential condition (I2.4bjl on M a(i is dM a p = (M aa C^ + M^C^du 1 , so that dM af3 = 
as a result of (|3,10bp . It is easy to check that (|3,10b|) implies (|3.10ap . It follows that there 
is a Lagrangian for (|3.9p with multiplier M a p if and only if M a p is constant and M a $ satisfies 
equation (|3.10b|) . Moreover, the Lagrangian is given by 

1 



L 



6 



M a/3 [3u^u 



2C* T v?u\ 



0- 



We remark that (|3. 10b|) is exactly the same as the condition for the existence of a bi-invariant 
symmetric bilinear form for a Lie algebra q with structure constants C" 7 . If q is semi-simple, 
the Killing form provides us with a non-degenerate solution to (|3.10bp . Consequently, (|3.9p is 
variational whenever q is semi-simple. Moreover, the number of solutions to (|3.10bh is equal to 
the dimension of the Lie algebra cohomology space #3(0). If 5 is simple, then dim #3(5) = 1 and 
the Killing form determines the only non-degenerate solution to (I3.10bp up to a scalar multiple 
(See [E], Theorems 11.1, 11.2). We remark that semi-simplicity is not a necessary condition for 
(I3.10bp to hold, as there are solvable Lie algebras which also admit bi-invariant bilinear forms. 
For example, consider the solvable 4-dimensional Lie algebra q of consisting of real matrices of 
the form 

CL\ (22 
03 CI4 


It is easy to check that the bilinear map M : q x q — > R defined by 

M(A, B) = Aa 3 & 3 + /2(a 2 &3 + 03^2 - aih - a^h] 
is non-degenerate and bi-invariant for all /i 7^ and all A £ R. 



4. Classification of Variational Systems in Two Dependent Variables 
In this section we establish a result characterizing the variational systems 
(4.1) u% y = f a (x, y, u 7 , u^, u^), 

of two equations and two dependent variables that admit multiple Lagrangians. In order to 
proceed, we need to make precise two concepts that are paramount to our discussion. We first 
define what it means for a system to have multiple Lagrangians. Then we review the notion of 
contact equivalence of two systems of differential equations (|4.ip . 

We say that a system of differential equations (|4.ip admits k Lagrangians if there exists a set of 
linearly independent Lagrangians {L\, L2, . . . , L/~} and a set of linearly independent variational 
multipliers {M^, M% p }, such that E a (Li) = M l a/3 (u^ y - p) for i = 1, . . . , k. We say that 
two /-Gordon systems ()4.ip are contact equivalent if there exists a local diffeomorphism 

$ : J 2 (R 2 ,R m ) -► J 2 (R 2 ,R m ), 
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where J 2 (R 2 , R m ) denotes the second-order jet-bundle of local sections s : R 2 — > R m , such that 

**(^-/ a ) = Q>2y-/ 7 )> 
and <I>*C C C, where C is the ideal generated by the 1-forms 

du a - u^dx — Uy dy, du" - u% x dx — u" y dy, du y - u" y dx - u yy dy. 

It was shown in |5j that any contact equivalence $ of two systems of the form (|4.1|) is the 
prolongation of a fiber preserving transformation 

(4.2) x = A{x), y = B(y), u a = C a (x, y, u"), 

up to an interchange x <-> y of the independent variables. 

We have the following theorem that completely characterizes the variational /-Gordon systems 
(|4.ip in two dependent variables which admit two or more inequivalent Lagrangians. 

Theorem 4.1. Let TZ denote the system 

(4.3) u xy = f(x,y,u,v,u x ,v x ,u y ,Vy), v xy = g(x,y,u,v,u x ,v x ,Uy,v y ). 

1) TZ admits three Lagrangians if and only if TZ is contact equivalent to a system 

u xy = A(x, y)u, v xy = A(x, y)v. 

2) TZ admits two Lagrangians if and only ifTZ is contact equivalent to a system 

u xy = W v (x,y,u,v), v xy = W u (x,y,u,v), 
where W satisfies one of W uu + W vv = 0, W uu = W vv , or W vv = 0. 

Remark 4.2. According to Proposition (|2.2j) and Remark (12. 3p . a system (|4.3p admits at most 
three Lagrangians. Moreover, since the dimension of the vector space of symmetric m x m 
matrices is m(m + l)/2, a system (|4.ip of m equations admits at most m{m + l)/2 Lagrangians. 
If a system of m equations (j4.1[) admits the maximal number of Lagrangians, then it can be 
shown the given system is contact equivalent to a system u" y = X(x,y)u a . 

The proof of the Theorem (|4.1j) depends on the following lemma that completely characterizes 
the /-Gordon equations u" y = g a (x, y, u 1 ) up to contact equivalence. The proof of Lemma (14.3|) 
is quite tedious and is delayed until after the proof of Theorem (14. lj) . 

Lemma 4.3. A system of partial differential equations u% y = f a (x,y,u 1 ,ui,Uy) is contact 
equivalent to a system u" y = g a (x,y,u y ) if and only if H = K and = 0. If the number of 
dependent variables m > 1 and H = K = AI, then A = X(x, y) and the given system is contact 
equivalent to a system u" y = \(x,y)u a . 

Proof of Theorem 1^.1. We first show that if a system (14. 3j) admits multiple Lagrangians, then 
(14.31) satisfies the hypotheses of Lemma (|4.3p and is contact equivalent to a system of the form 

(4.4) u xy = F(x, y, u, v), v xy = G(x, y, u, v). 

According to Proposition (12. 2p . if (|4.3p admits a first-order Lagrangian L with a symmetric 
multiplier M a p, then M a/3 satisfies the algebraic conditions 

(4.5) M^HJ = MpylQ, M aa S^ + M 0a S^ = 0, 

where ilff, K2, and are given by (|2.3p . Multiplying the second equation of (|4.5|) by M a/3 
yields 

(4.6) = M af) M aa S° + M af3 M Pa S° = 5^ + ftS* = 2S" 
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If their are only two dependent variables, then equation (|4.6p implies that 



qa _ 

expressed as 



(4.7) 



12 - "12 - °- SinCe 

—S® g , we have SfL = for all a, (3, 7 = 1,2. Then the algebraic conditions (|4.5|) can be 



02 

°12 



Mi 2 

M 22 



0. 



Hi - Kl 


Hi 


-Kl 








Hi 


-Kl 


Hi - Kl 


Hi 


Hi 


-Kl 


-Kl 


Kl 


Kl 


-Hi 


-Hi 



If (|4.ip admits 2 or more inequivalent Lagrangians, then there are least two linearly independent 
solutions to (|4.7p and it follows that the rank of A is at most one. Moreover, the rank of A is 
one or less if and only if H = K and is rank zero if and only if H = K = XI. From Lemma (|4.3|) . 
we deduce that Rank A < 1 if and only if (|4,3p is contact equivalent to (|4.4p . 

We will complete the proof of the theorem by analyzing the algebraic conditions (|4,7p for 
systems of the form (|4.4p . A calculation of H and K for (|4.4p reveals that 



H = K 



' F F 

G u G v 

In this case the algebraic conditions (I4.7P for the existence of a multiplier M af3 simplify to 

(4.8) F V M U + (G v - F U )M 12 - G U M 22 = 

The differential condition (|2.4bp reduces to dM a p = 0. Consequently, solving the variational 
multiplier problem for (|4.4p is equivalent to determining all constant solutions M a $ to the 
equation (|4.8p . 

There are 3 linearly independent solutions to (|4.8p if and only if H = K = AI. In this case 
(|4.4p is contact equivalent to a system 

(4.9) u xy = X(x,y)u, v xy = X(x,y)v. 

The most general Lagrangian for (|4.9p is a linear combination of the Lagrangians L\ = u x u y + 
Xu 2 , L 2 = u x v y + Xuv, and L3 = v x v y + Xv 2 . 

We analyze the case where the rank of A is exactly one and we assume there are two non- 
degenerate, linearly independent, constant solutions Mig and M 2 ^ to (|4.8p . We claim there 
exists an indefinite multiplier M = (M aj3 ) that satisfies (|4.8p . If one of det Mi < or detM 2 < 0, 
then we are done, so we assume that det Mi > and det M 2 > 0. If 

b 



(4.10) 



Mi 



M 2 



then it follows that ac > and pr > 0. We claim there is a scalar /x such that det (Mi - 
Indeed, if 



-uM 2 ) < 0. 



(detM 2 ) 2 u 2 



+ (2bq 



(4.11) P(u) = det(Mi - «M 2 ) 

then the discriminant A can be expressed as 

A= {aq- bp) 4 + 2(p 2 d et Mi + a 2 det M 2 )(qg 



ar 



pc) u + (det Mi 



6p) 2 + (a 2 M 2 -^ 2 Mi) 



a 2 p 2 

It follows from (^TUj) and (|4TT2|) that A > with equality holding if and only if Mi = t M 2 for 
some t € R. Consequently, the polynomial (|4.1ip has two real roots and there exists /x E R such 
that det (Mi — «M 2 ) < 0. Now we have established that if there are two independent solutions 
to (|4.8p . with at least one of the solutions non-degenerate, then there is an indefinite multiplier 
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If we make a linear change of variables u a = T^u 1 ', then a direct calculation of the Euler- 
Lagrange equations for the Lagrangian L and the transformed Lagrangian L verifies that the 
corresponding variational multipliers transform according to the rule 

(4.13) M ap = T°Tp M aT . 

After a linear change of variables u a — > T"ii 7 , we may then assume that the indefinite multiplier 
M a p is of the form 



(4.14) M 



1 

1 



Substituting (|4.14p into f|4.8j) implies that G v — F u = 0. As a consequence of the de Rham 
theorem, we see there exists a smooth function W(x, y, u, v) such that 

(4.15) W v = F, W u = G. 

There is a second multiplier N, independent of M, satisfying (|4,8p . We may assume, possibly 
after subtracting a scalar multiple of M, that 

N=( a ° 
\0 b 

According to (|4.13l) . a linear transformation u — > Xu, v — > A _1 v preserves M and transforms N 
as 

'aA 2 \ 
6A- 2 j ' 



N 



We may then assume, possibly after a scaling N — > fcN, that N has the form 

T N 



where £ = 0, 1, or —1. Taking into account that (|4.15p holds, substituting N into (|4.8p implies 
that 

W vv - eW uu = 0. 

If £ = —1, then W satisfies Laplace's equation W uu + W vv = and there exists a function 
Z(x,y,u,v) such that Z u = W„ = F and — Z w = W u = G. The most general Lagrangian in this 
case is given by 

L = Cl(U X Vy + W) + C 2 {U X Uy - V X Vy + Z) . 

If e = 0, then W OT = and (|4.4p can be expressed 

(4.16) u xy = a'(u), v xy = a"(u)v + b'(u). 

The most general Lagrangian for (|4.16p is 

L = c\\u x u y + 2a(u)] + C2[w x «i/ + a '( u ) v + ^( u )]- 

If e = 1, then W satisfies the wave equation and W = W\(u + v) + W2(u — v). The most general 
Lagrangian in the case is 

L = c\[u x u y + v x v y + 2W\{u + v) — 2W 2 (« — v)] + C2[u x v y + W\(u + v) + W2{u — v)\. 

This establishes the second statement of the theorem. □ 
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Proof of Lemma 4-3- If we assume that two systems u xy = f a and u xy = f a are contact 
equivalent with the change of coordinates given by (|4.2|) , then the formulas ()2.3f) and a tedious 
application of the chain rule will verify that the transformation rules for H, K, and S are 

(A-17) H a =^— — —H a K a =^— — —K (T g« = d ^ dyT du " o« 
1 ' * A'B' dvF du@ ' A' B' du a duP ' ^ dvP dvP dvH Ta 

where (du a /du 1 ) ■ (du 1 /du 13 ) = 8". It follows from (|4.17p that the conditions H = K and S = 
are invariant with respect to transformation (j-4. 2f) . 

For a system u xy = g a (x,y,u' y ), using (12.3f) we see that 

(4.18) H " = K "= d £i> S ^ = °- 

From (|4.17p and (|4.18p . we deduce that any system u xy = f a (x,y,u y ,u2,u y ) that is contact 
equivalent to u xy = g a (x,y,u y ) must necessarily have H = K and S^L = 0. 
We now prove that any system 

(4.19) u« y = r(x,y,u\ul,ur) 

with the property that H = K and S9[ = is equivalent to a system u x = g a (x,y,u y ). It 



follows from flU} that if H° =K™, then 



d 2 f a d 2 f a 

°j „ ^ r = °> Va,/3,7 = 1,2, ... ,m. 



duxdux duydu'y 



Consequently, equation (|4.19|) simplifies to an equation of the form 
(4.20) u% + C a pi vP x ul + + B°u] + G a = 0, 

where C" 7 , ^4", and G a are functions of and u e . Moreover, S 1 ^ = if and only if 
C" 7 = C~jg. With a judicious choice of coordinates, we will now eliminate the quadratic terms 
of (|4.2U|) . If we let u a = g a (x,y,u" / ), then (|4.2U|) transforms as 

(421) w< + - c - {x ^ gT) ^p%) + A >* + 5 > l + Qa - 

We see that CfL = whenever g a satisfies 

< 4 - 22 > ^ +c -^-" T >^% = 

We differentiate (|4.22|> with respect to u* 5 , and after substituting from (14.221) and skew-symmetrizing 
over (3 and 8, we obtain the integrability conditions on (|4.22p 

On the other hand a calculation of H and K for (|4.20p yields 

B 2 BC a BC a 

(4 ' 24 > 8^ iH " - K °1 = 8# - ft? + - C ^ 

As a consequence of (|4.24p . the integrability conditions (|4.23p are satisfied whenever H = K. 
The system of partial differential equations (|4.22p then satisfies the Frobenius condition and we 
deduce that there exists, at least locally, a non-degenerate collection of functions g a satisfying 
(Q2| . 
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We may now assume that u" y = f a {x,y,u 1 ,ul,Uy) is of the form 
(4.25) u% + A«(x, y, u e )u2 + B%{x, y, u> 7 + G a (x, y, u e ) = 0. 

For (|4.25p . we calculate 

dA° . dBfi „ dA% dB* 

("•) «? - *? = ^ + ^t«S + W - ^ + ^ - -s^- 

Evidently, iJ^ = K% only if dB^/du = dA^/du 13 = for all a, (3, and 7. It follows that (|425j) 
simplifies to 

(4.27) u% + A»(x, y)n2 + B«(x, y)u} + G a (x, y, u") = 0. 

We now eliminate the functions A"(x,y) and B"(x,y) from (|4.27p with a transformation u a — > 
iV!j*(x, y)u 7 , where satisfies the system of partial differential equations 

8N a 8N a 

(4.28) _JL + ^ = 0l ^L + B?N; = 0. 

The integrability conditions for (|4.28p are given by 

(4.29) _L + A r B c* = _L + B T A a_ 

As a consequence of (|4.26p we see that (|4.29p holds whenever H = K. It follows that the system 
(|4.28p satisfies the Frobenius condition and there exists functions N"(x, y) such that det 7^ 
and (|4.28|) is satisfied. Moreover, we arrive at an equation of the desired form 

(4.30) u% = g a {x,y,v?). 

To prove the final statement of the lemma, we assume that (|4,30p has the property that 
= K" = X5°. A routine calculation shows that H" = K° = dg a /du 1 . It is evident that 
= dg a /du 1 = X5" if and only if dg a /du 13 = whenever f3 ^ a. This means that each 
function g a satisfies 

(4.31) ^ = A, a = l,2,...,m, 

with no summation on a. Differentiating equation (|4.3ip with respect to u 13 , 1 < (3 7^ a < m, 
results in dX/du 13 = for all (3. We deduce that A = X(x,y) and it becomes apparent from 
(I4.3ip that g a = X(x,y)u a + k a (x,y). After a transformation u a — > u a + m a (x,y), where 
m xy = ^ mCX + we see that (|4.30p is equivalent to u% y = X(x,y)u a . □ 

5. The Variational Bicomplex For Systems of PDE 

In this section we introduce some basic definitions and results used in our solution to the 
variational multiplier problem in Section 6, including infinite jet bundles and variational bi- 
complexes. As we are only interested in applications to our study of the variational multiplier 
problem, our discussion will be of a rather brief nature. For a detailed and intrinsic construction 
of variational bicomplex, we refer the reader to p], [3], [18], and |25j . 

Let TT k : J k (E) — > R n denote the bundle of fc-jets of local sections of the trivial bundle 
E = R n x R m . Local coordinates for J k (R n ,R m ) are given by 

where 1 < i\ < 12 < • • • ik < n and 1 < a < m. There are natural projections ni : J l (E) — > J k (E) 
for / > k. The infinite jet bundle over E, J°°(E), is defined as the inverse limit of the sequence 
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d H f = (A/) dx\ d v f = -zL.fr + ^e? + e« + • 



of finite jet bundles {J k (E) \ k = 0, 1, 2, . . . }, along with the projections 7r£° : J°°(E) — > J fc (£) 
and vr 00 : J°°(i?) -> R n . The contact idea? C(J°°(E)) is generated by the 1-forms 

Ofxiz-ik = ^ U hi2---ik ~ U hi2---ikj^ xJ ' V = 0, 1, 2, . . . . 

The full exterior algebra f2*(J°°(.E)) of differential forms on J°°(E) is generated by the 1-forms 

dx\ # a , Of, dfj, ... . 
There is a bi-grading of the differential forms on J°°(E), 

n p (j°°(E))= n r ' s (j°°(E)), 

r+s=p 

where £l r ' s ( J co {E)) is the C°°( J°°(£'))-module generated by differential forms of the type 

dx 11 A dx h A • • • A dx ir A 0% A • • • A 0£. 

The exterior derivative d : n p (J°°(E)) — > fF +1 (J°° splits into the horizontal and vertical 
differentials d = du + dy , where 

d H : n r ' s (J°°(E)) fi r+1 ' s (J 00 (£;)), : n r ' s (J°°(E)) fi r ' s+1 (J°°(i?)). 

Since d 2 = 0, it follows that d 2 H = 0, d v = 0, and djjdy = —dydn- The local coordinate 
expressions for the horizontal and vertical derivatives of a smooth functions / £ C 00 (J 00 (E)) 
and 1-forms dx 1 and 0" are given by 

(5.1) d#0? = dx 1 A 0£, d y 0? = 0, dnidx 1 ) = 0, dyidx*) = 0, 

where A denotes total differentiation with respect to x l . The free variational bicomplex is defined 
to be the double complex {Q r,s J°°(E), dfj, dy } s >0; r=o,i,...,n- 

In our solution to the variational multiplier problem for the /-Gordon systems 

(5.2) u% y = f a (x,y,u f ,ul,u' y '), a = 1, ...,m, 

we investigate the existence of certain cohomology classes in the constrained variational bi- 
complex associated to a system of partial differential equations. To construct the constrained 
variational bicomplex associated to (15. 2|) , we begin with a trivial bundle tt : R 2 x R m — > R 2 and 
consider the second-order jet bundle J 2 (R 2 ,R m ) with coordinates given by 

(x, y, u a ,u*,u*,u* x ,u* y , u% y ), a = l,...,m. 

An /-Gordon system ()5.2[) defines a (5m + 2)-dimensional submanifold 1Z 2 A J 2 (E) called the 
equation manifold of (|5.2[) . We define the first prolongation of 1Z 2 as the 7m + 2-dimensional 
submanifold TZ 3 A J 3 (£) defined by fl52J and = D x f a and = D y f a . Further dif- 
ferentiations of (15. 2D will yield submanifolds 7£ fc — > J k (E). For convenience we define 7?.° = 
and TZ 1 = J 1 (S). We define the infinite prolonged equation manifold 7Z°° to be the inverse limit 
of the sequence {7Z k \k = 0, 1, 2, ... }, along with the natural projections : 7Z°° — > M and 
7r^° : 7£°° — ► 7£ fc . We remark that there is a unique map : 7Z°° — > J°°(E) that satisfies the 
commutative diagram 

7e°° — J°°(E) 



1Z k — ^ J k (E). 
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For an /-Gordon system (j5.2[) . coordinates on 7Z°° are given by 

/ GLOLCLOLGLQ.Q.GL OL \ 

\X,y,U , U x , Uy , U xx , Uyy j V'xxx i ^yyy > ^xxxx i ^"yyyy >•••/• 

We define the contact ideal C(TZ°°) on 1Z°° via the pullback of the contact ideal on J°°(E), that 
is C(TZ°°) = l*C(J°°(E)). The contact ideal C(ft°°) for an /-Gordon system ([5J2D is generated 
by the 1-forms {9, 9%, 9%, 9 xx ,9^ y ,9^ xx , 9° yy , ...}, where 

du a — u x dx — u y dy, 9 X = du x — u xx dx — f a dy, 9 y = du y — f a dx — u yy dy, 



&xx = du*x - Kxxdx - i*(D x f a )dy, 9^ y = du^ y - i*(D y f a )dx 

A basis for the C 00 (7?. oc ')-module of 1-forms r2 1 (7^°°) is given by 

(5.3) {dx,dy,9 , 9 X , 9 y , 9 XX , 9 yy , 9 XXX , 9 yyy , . . . } 

For every p, we have a bi-grading S7 P (7£°°) = Q r ' s (TZ ao ), where r + s = p and J7 r,s (7?, 00 ) = 
i^fi r ' s (J°°(^)). The exterior derivative d : W(R,°°) -> ftP +1 (ft°°) splits as d = d H + dy, 
where dn '■ Q r ' s — > f2 r+ ' s and dy '■ Q r,s — ► O r,s+1 are the horizontal and vertical differentials, 
respectively. In the following section we will frequently use the djj and dy structure equations 
for the coframe (|5.3|) . which are given by du dx = du dy = dy dx = dy dy = 0, 

(5.4) 



"XX 

yy 



ytjii 



dy,... 



d H 9 a = dx A 0° + dy A 0£, d H 0" = dx A 9% x + dy A dy f 
d H 9 a y = dx A dyf a + dy A 9% y , 
d H g = D x gdx + D y gdy, dy g 



d v 9 c 
d 9 no 



dy9 xt 



d v B c . 



y J 



0, 



+ 



du c . 



where g € C 00 ^ 00 ) and D x and -D y denote total differentiation constrained to the equation 
manifold 7£°° . 

Definition 5.1. The constrained variational bicomplex {l*'*(TZ°°, dn , dy) associated to an /- 
Gordon system is the pullback of the free variational bicomplex Q*'*(J°°(E),dH,dy) by too : 







R 



^(ft 00 

civ 



ft°' (ft° 



dy 



d n 



dv 



d v 

SI 2 - 2 (ft 00 ) 

dy 
dy 



The columns of the variational bicomplex on 1Z°° are locally exact, while the rows will not be 
exact in general. We define the horizontal cohomology classes of the variational bicomplex by 

{u g n r ' s (n°°)\d H uj = o} 



H r > s CR°°) 



1,2, s>0. 



{d H v\v € n r-1 '*(W°)} ' 
It is easy to see that H r ' s (TZ°°) is a vector space over R. The cohomology classes in H r ' s (lZ°°) 
have interesting interpretations. For example, each class [u] = H 1,0 (TZ°°) represents a classical 
conservation law. Indeed, if oj = M dx + N dy, then dfiw = if and only if D X M = D y N when 
restricted to the equation manifold 1Z°°. We refer to a cohomology class [to] £ H l,s {JZ°°) as a 
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type (l,s) conservation law or form-valued conservation law. In the following section, we will 
show that the solution to the variational multiplier problem is closely related to the existence of 
non- trivial classes [u] G H^ 2 (n°°). Since systems of the form ()5.2[) are of Cauchy-Kovaleskaya 
type, it follows from a general result of Vinogradov [25] that the horizontal cohomology spaces 
H Q < a (K°°) satisfy 

(5.5) H^CR 00 ) = R, H^CR 00 ) = 0, s > 0. 

We remark that (15. 5ft was also established in [5j by constructing a coframe adapted to systems 
(|5T2l) . 



6. Derivation of Necessary and Sufficient Conditions for the Existence of a 

Variational Multiplier 

Our first result states that the problem of determining all Lagrangians and variational mul- 
tipliers for an /-Gordon system is equivalent to determining all d closed forms u G Q 1,2 (7?. 00 ) 
of a certain type. We also give a description of the general form of possible Lagrangians for 
an /-Gordon system. Finally, we show that a variational multiplier has no one-jet dependence. 
Proposition (|6.ip along with Corollaries (16. 2p and (16. 3p will suffice to establish Proposition <|2.1|) . 
which was stated without proof in Section 2. 

Theorem 6.1. For a system of differential equations uf. y = f a (x,y,u r ,ul,Uy) the following 
statements are equivalent. 

(i) There exists a type (1,2) form 

(6.1) oj = {T a0 dx + S a/3 dy) A a A 6 P + R a0 {6 a A 0% A dx - 0? A 0£ A dy ) 

such that du = on TZ°°. 

(ii) There exists a first-order multiplier M Q/3 (x, y, u 1 ,u2, viy) and a first-order Lagrangian L(x, y, u 1 , u2, Uy) 
such that E a (L) = M af3 (u^. y - f 3 ). 

(Hi) There exists a multiplier M af3 = M a p(x,y,u r ) and a Lagrangian 

L = -R a0 (x, y, u 7 )u"n^ + Q a (x, y, u 7 )u" + P a (x, y, u 1 )^ + N(x, y, u 1 ) 
such that E a (L) = M afj [uP xy - f ) . 

Proof. We first show that (i) implies (iii). Suppose that uj is given by (|6.ip and that duj = 0. on 
TZ°°. It follows immediately that djjoj = dyuj = 0. A routine calculation using (15. 4|) shows that 
if dyuo = 0, then R a0 = R a p(x, y, u e ) and u> is of the form 



1 fdR, 



OR, 



+ 



dR^p 



du c 



^ul + T^y^) 



du° 



^f) Uy< + S o a0 (x,y,u e ) 



9 a A 9 A dx + R a0 9 a A0 f3 Adx 
6 a A 6 13 A dy - R pa 6 a A 0£ A dy. 



If we define Po G fl^^Tl 00 ) by 

Po 

then using (|5.4p we see that 



-Rapul 9 a A dx + Rp a vP y 9 a A dy, 



to - d vPo = t° e a a e fJ a dx + SL e a a , 



A dx. 
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Since dyiuj — dypo) = 0, there exists a form p\ € S^ 1 ' 1 (72.°°) such that dyp\ = uj — dyPo- Since 
the functions S^p and have no 1-jet dependence, we may choose p\ to be of the form 

pi = -P a {x, y, u e ) 6 a Adx + Q a (x, y, u £ ) 6 a A dy. 

We now define p = po + p\ so that dyp = uj. Since dydji = —dydn, we have dy(dup) = 
—d}jdyp = —dfjuj = on 1Z°°. Therefore there is a Lagrangian form A £ il°' 2 (7^ 00 ) with 
dy\ = diip on TZ°°. We will now show that A = Ldx A dy, where L is of the form given in 
statement (hi). Computing dnp on J°°(E) yields 

d-HP = ~ [{Rafi + Rpa) u p xy + g a (x, y, u e ,u x ,u e y )] 6 a Adx A dy- 
(Rapu? + P a ) 9y Adx Ady — (R pa u^ + Q a ) 6 X Adx A dy, 

where 

'dRp a dP a \ dP a dQ a 



dx du 13 J y dy dx 

When restricted to 1Z°°, 

-d H p = [(2R {a0) f + g a ) 9 a + (R a pu r i + P a ) 6° + (R 0a u^ y + Q a ) 0*\ Adx A dy, 

where Ri a p) = {Rap + R/3a)/2 and g a is given by (|6.2p . Since A = Ldx Ady and dy\ = dnp on 
1Z°°, we see that L must satisfy 

dL dL dL 

— = -2R (a0) f -g a , — = -R^vPy -Q a , — = -R a0 u p x - P a . 

It follows that there exists a function N(x,y,u e ) such that 
(6.3) L = — {Rapu 13 ^ + Q a u a x + P a v% + N(x, y, u e )) . 

If we apply the Euler-Lagrange operator E(\) = E a (L) 9 a Adx A dy on J°°(E), then 

d H p + E(X) = d v (X). 

Since dnp = dyX when restricted to the equation manifold TZ°°, we deduce that implies l*E(X) = 



0. On the other hand, a direct computation of the Euler-Lagrange equations for (|6.3p gives us 

dL 

E a (L) = 2R {a/3) u x \ y + 9a + Q^- 

Since i*E a {L) = for all a, we get g a + dL/du a = — 2iZ( Q(3 )/ /3 , which implies 

E a (L) = 2R {a3) (4 y -f), 

and (hi) is proved. Moreover, the multiplier M af3 is defined by M af3 = R a/3 + Rp a . 

We now prove that statement (hi) implies (i). Assume that there is a variational multiplier 
M a p(x, y, u 1 ), a hrst-order Lagrangian X = Ldx A dy, where L is of the form (16. 3h . and 

(6.4) E a (L) = M af) (ul y -f). 

Note that (|6.4p explicitly determines that the multiplier is M a p = R a/3 + Rp a . By the first 
variational formula (See [I], Corollary 5.3), we have E(X) + dfjT] = dyX, where 

r)T f)T 
7] = — 9 a Adx — — 9 a A dy. 
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Since l*E(X) = 0, it follows that t*dnr] = dyX. Then define u = dyrj and the resulting calculation 
of du on the equation manifold is 

du) = dffdyOJ + dyLU = dfj{dyr]) + dydyr\ = —dydyX = 0. 

Moreover, using (|5.4p to compute dyq will verify that uj is of the form (|6.ip . 

We have proven that (i) is equivalent to (iii), and clearly (iii) implies (ii). Assume that 
(ii) holds and there is a first-order Lagrangian L(x,y,u' y ,u2,Uy) and a variational multiplier 
M a p(x, y, u J ,u2, Uy) for the system u^ y = f 13 . We will show that this implies (iii). Calculating 
the Euler-Lagrange equations for an arbitrary first-order Lagrangian, we find that 

ran vfT\- d 2 L p ( d 2 L d 2 L \ p d 2 L 

(6 - 5) Ea{L) ~ ~o^d7 x Uxx " [d^d4 + duJ^E ) Uxy " d^? y Uyy + Ga ' 

where G a is a first-order function. On the other hand, we are assuming 
(6.6) E a (L) = M ap « - f(x, y, u\ u%,u$)) , 

where M a p is a non-degenerate first-order multiplier. Comparing (|6.5|) and (|6.6|) . we deduce 
that 

(0.7) J^ = 0, =a 

It follows that L is of the form given in statement (iii) and the fact that M a p has no one- 
jet dependence follows immediately from a calculation of the Euler-Lagrange equations for a 
Lagrangian satisfying (|6.7p . □ 

The following corollary states that a Lagrangian A corresponding to a type (2,1) cohomology 
class [uj] is unique up to modification by a total divergence. 



Corollary 6.2. Let uj be given by (16. ip with du = 0. If 

X = L(x, y, v? ,uj, u^) dx A dy, A' = L'(x, y, u 1 , uj, it^) da; A dy 
are two Lagrangians and 77, 7/ G J7 1,1 (7?. 00 ) satisfy 

dnV = dy\, dHT)' = dyX', dyn = dyn' = uj, 
then A' = A + d H j3 for some form (3 G Q lfi {TZ°°). 

Proof. Since dy(n' — 77) = 0, there exists a form /?o G fi 1 ' ^ 00 ) such that rj = 77 + dyfio. It 
follows that 

d^A' = dnrf = dniv + dy/3 ) = dyX + dudyPQ = d v (X - dnPo)- 
Consequently, A and A' satisfy 

d v {X' -X + d H (3 ) = 0, 

or equivalently that 

A' — A + dnflo = a{x, y) dx A dy. 
Defining j3 = /?o — A(x, y) dy, where A x (x, y) = a(x, y), we have that 

77 = 77 + (fy/3 and A' = A — dn/3, 
as required. □ 

The following corollary gives a description the general form of a /-Gordon systems that is vari- 
ational. Corollary (|6.3p . along with Theorem (|6.ip and Corollary (|6.2p . establishes Proposition 
(O). 
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Corollary 6.3. If u@ = f^(x, y, u 1 , u2, Uy) has a nonsingular variational multiplier M a p, then 

f = Ciu%ul + A*u% + Blu\ + CP, 
where A^,B^, Cf r , G G C°°(J°(E)). 



Proof. According to Theorem (|6.ip , we may assume that there is Lagrangian L of the form (|6.3p 
and a variational multiplier M a p(x,y,u), with detM a/3 ^ and 

(6.8) E a (L) = M ap ( U P y -f). 

On the other hand, the Euler-Lagrange equations for (|6.3|) are explicitly given by 

(6.9) E a (L) = (R a , + R, a )ul y + + ^ " ^) <ul + 

( dQg dQ p 8R a A g (dP a dP a dR a0 \ dN dQ a 8P a 
V dvP du a dy J Ux \ duP du? dx ) Uy du a dx dy ' 

Clearly, M ap = (R a p+Rp a ) and multiplying equations ([62]) and (E9j) by M" 7 , where M° 7 M 7/3 = 
5p, yields the desired result. □ 



In view of Theorem (16, ID , we see that solving the multiplier problem is equivalent to finding 
all type (1,2) forms u of the form (16.11) with dw = 0. From Corollary (16. 3|) . we can restrict 
ourselves to examining systems of partial differential equations 

(6.10) u% + C« e (x, y, u^ulul + A«(x, y, u p )ul + B%(x, y, u> 7 + G a {x, y, v?) = 0. 

For systems (|6.10p . the following proposition gives necessary and sufficient conditions for dw = 0. 

Proposition 6.4. There exists a differential form uj G il 1 ' 2 (7?. 00 ) of type (I6.ip with duj = if 
and only if M a p = (R a p + Rp a ) /2 satisfies M a p = M a/3 (x,y,u y ) and 

(6.11) M ai m = M 01 Kl, 



3-12) M ai Sl = -M 01 S2„ 



l al Up — lVlp-f . 

(6.13) dM a0 = M ai Vfp + M Pl Q.% , 
where Sj t = Cj e — and Qa = Cj e du € + dx + Aa dy. 

Proof. In the proof of Theorem (|6.ip we showed that if dyoj = 0, then R a p has no 1-jet depen- 
dence. Since M a/3 = (R a p + R/3 a )/2, it follows immediately that the functions M a p depends 
only x, y and u. 

To simplify our calculations, we define uJ = oj — jdn {Rap — Rpa) 9 a A9 13 . A routine calculation 
using (15. 4p yields 

(6.14) J = {T ap dx + V ap dy) A 9 a A B & + M a0 {9 a A 9% A dx - 9? A 9% A dy), 

where M aj3 = {R a p + R /3a )/2. Without loss of generality we may assume that T a p and V a p are 
skew-symmetric. Using (15. 5p and the exactness of the columns of the variational bicomplex, it 
can be shown that dw = if and only if dnuj' = and dyw' is da exact. We will later show 
that if djjw' = for (|6.14p . then dyco' is necessarily dfj exact. 

We first calculate djiuj' and show that duuJ 1 = if and only if (I6.1ip - (l6.13p hold. Using (15. 4j) 
we calculate 

(6.15) d a J = [(D x V a p - D y T a0 ) 9 a f\9 p - 2M a0 9 a A (d v f)] AdxAdy- 

[(D y M a p - 2V a0 ) 9 a A9? + (D x M a p + 2T a0 ) 9 a A 0$\ A dx A dy. 
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It follows that djfUj' = if and only if 



(6.16) 



(6.17) 



D x M a0 + 2T Q0 + 2M, 



dp 



a/3 



D x V a/ 3 — D y T a0 



0, D y M a(3 -2V a(3 + 2M aj ^j = 0, 



<9u c 



0. 



Decomposing the two equations (|6.16p into their symmetric and skew symmetric parts produces 
four equations 



(6.18) 
(6.19) 



D m M. 



dP 



dP 



L a/3 



M, 



01 



a/7 

du% 



M, 



y 

dP 



DyM af3 



-M, 



M 



ai d4J " 2 



dp' 

dul "*d& 
dp dp 

M ^d7 x ~ M ^ 



Substituting (|6.19p into (16.171) and taking (j6. 18f) into account, we arrive at 

(6.20) M ai (Hl + Kj) = M^HZ + 

where 



W 

- L - L r\ 



dp , df 



+ 



du a du T y du% 



Dr 



dP_ 

du* 



dp dp 



du a du T x du* 



dp 
dvS, 



It follows that dftio' = if and only if (|6.18p and (|6.20p hold. For a system of differential 
equations (|6.10p . we expand (|6.18p to arrive at 



BM, 



ap dM a0 



dx 



+ 



du" 



4 = M a7 [C^ul + Bj) + M M {Bj a u% + Al) , 



dM a0 dM a0 



Oy 



+ 



du- Uy 



M a , (Cj e ul + AJ) + M M [Cl e u% + Al) . 



Since M a0 has no first-order dependence, we have 

dM a0 _ ^ ^ , ^ ^ dM af3 



(6.21) 



0M e 

dx 



M ai Bj + M 0y B2 



du e 
dM a0 

dy 



M a ^ + MM 



M^Al+MfrAZ. 



It follows that (|6.2ip is equivalent to conditions (|6.12p and (|6.13p . If (|6.12p holds, then it can 
be shown that the integr ability conditions d 2 M a p = for (|6.13p are 



(6.22) 



M aa {H% - Kf) + M Pa {HI - Kl) = 0. 



Equations ([6T22]) and (f6T2Q]) both hold if and only if M aa H% = M 0a K°, which is precisely 
condition (|6.1ip . 

We now show that dyu' is dn exact whenever djiu 1 = 0. We claim that dn( = dyw' for 
(6.23) ( = ±M a<7 S^e a AO? A9~< 

From (|6.10p . (|6.14p . and (|6.19p . we deduce that u' can be expressed as 

J = M aa (C° u T x + Bp) 6 a A 9 P A dx + M a0 Q a f\Q x l\dx — 

M aa (c° T u T y + a%) e a a e* 3 a d y - M a0 e a a el a d y . 



THE INVERSE PROBLEM OF THE CALCULUS OF VARIATIONS 21 

We now calculate dyu>' — djj^. Since we are assuming that dnio' = 0, we use (j6. 12|) . (|6.13p . and 
the fact that = C^ 7 — to arrive at 

(6.24) d v J - d H ( = Q aPl 9 a A 6 P A # 7 A dx + R aPl 9 a A 9 P A # 7 A dy + 
[{\M aa S a 01 - M^C° afj ) 62 A dx - {\M aa S^ - MyrCZp) 9~< A dy] A 9 a A # /3 , 

where Qa/3j and i? a/ 3 7 are totally skew-symmetrized expressions depending on M a/ 3, CfL, A", 
and their derivatives. If we skew-symmetrize over a and /3, it follows from (|6.12p that (|6.24[) 
simplifies to 

(6.25) d v J - d H ( = Q a ^0 a A fi A # 7 A dx + i?^ 7 ^ Q A9 P AO 1 Ady 

Since dn{dyuj' — dnC) = — (dfl-u/) = 0, applying d# to the right hand side of (|6.25p produces 
the equation 

(6.26) [SQa^y - + {DyQafl-y - D x R a/3l )9 a ] A 9 P A 9~< A dx A dy = 0. 

It follows form (|6.26j) that Q a p^ = R a p<y = and we can conclude that dyu' = dnC- □ 



7. Concluding Remarks 

We conjecture that the methods used in this paper can be used to obtain a complete solution 
of the inverse problem for systems of the form 

(7.1) g l \x,y,v?)h a ^ + f a (x\u\ul) = 0, 

where a, {3 = 1, . . . , m and i, j = 1, . . . ,n. In particular, if det ^ and det h a p ^ 0, then 
the system (17.ip is a generalization of harmonic map equation. A partial solution of the inverse 
problem harmonic map equation was obtained was obtained by Henneaux |14| . It would be 
some interest to derive a set of invariant analogous to the invariant derived in this paper that 
completely characterize the existence of Lagrangians for systems (17. ip . 

Appendix 

To complete our solution to the inverse problem, we have the following general result on 
existence of solutions to systems of total differential equations with an algebraic constraint. 
Proposition (jA.ip establishes (|2.9p and (|2.1ip and completes the proof of Theorem (|2.4p . 

Proposition A.l. Consider the system of total differential equations 

dz a 

(A.l) — = A* d (x)z 1 1 < a < m; l<i<n, 

ox 1 ' 

coupled with the algebraic condition 

(A.2) B^{x)z^ = 0, 1 < a < I. 

Suppose the integrability conditions for (jA.ip are satisfied whenever (|A.2p holds and that given 
xq £ R n , there exists a neighborhood U of xq such that Rank{B® (xo)} = Rank{B!^ (x)} = k < I, 
for all x € U. We also assume that the algebraic system ()A.2p is complete in the sense that 
differentiating (|A.2p and substituting from (jA.ip produces no new algebraic conditions. Given 
xq E U C R n and zq 6 R m with B^(xq)zq = 0, there exists unique smooth functions z a = z a (x) 
defined on a neighborhood V C U of xq, such that z a (x) satisfy both (jA.ip and (jA.2|) . and 
z a (x Q ) = z$. 
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(A.3) -^Jr + B% A°i = G% B^. 



Remark A. 2. For the more geometrically inclined reader we note that Proposition (|A.lj) is 
a special case of the following result: For an m + n dimensional manifold M let X C Q 1 (M) 
denote an exterior differential system of rank m and let q be a regular value of a smooth map 
F : M -> Q. Assume for all p G S = F~ 1 (q) we have C Ker (F« : T p M -> T g Q) and that 
X. s is a Frobenius system. Then for all po G X C M there is a unique maximal n-dimensional 
integral manifold <f> : N —* M through po such that 0(-/V) C S. 

Proof. Since we are only constructing local solutions to (jA.ip and (|A.2p . we will assume that 
Rank{S^(x)} = k for all x G R n . Assume that no additional algebraic constraints are created if 
we differentiate ()A.2j) with respect to x 1 and substitute from (1A.1I) . This guarantees the existence 
of functions G%{x) such that 

dB« 

dx l 

The system of differential equations (jA.ip are associated with the C°° distribution A on 
R n x R m generated by the vector fields 

( A -4) ^ = ^ + ^ Z<T ^' i = 1 . 2 .-»"- 

We then define the C°° function F : R n x R m -> R' by F(x, z) = B*(x)z~< and let £ = F -1 (0). 
As a consequence of (|A.3[) we see that Rankj-Di^ ^} = Rank{B"(x)} = A; for all (x,z) G S. 
It follows immediately that S C R n x R m is regular submanifold of dimension m + n — k. In 
addition, we can use (|A.3P to show that X^B^z 1 ) vanishes identically for all (x,z) G S. We 
conclude that for any p G £ and i = 1,2, ... ,m, 

(A.5) G Ker(F, : T p R n+m ^ TqR') = T p £. 

It follows from (fOj) that Ai s is a C°° distribution on the submanifold S. (See Boothby [6], 
Lemma IV. 2. 4.) If the integrability conditions are satisfied for f| A. 1 1) whenever (|A.2p holds, then 
A|s is involutive. For any point (xq, Zq) G E, we can invoke the Frobenius Theorem to obtain a 
unique maximal integral manifold <p : V — > S with V C R n , 0(0) = (xo, zq), and T p 0(V) = A p . 
Moreover, <j) is C°° viewed as map (p : V — > R n+m and (j)(V) defines the graph of a smooth 
solution z a (x) to the system of differential equations (jA.ip . Since <j)(V) C S, we see that z a (x) 
satisfies the algebraic condition (|A,2p . □ 
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